Muon decay in a linearly polarized laser field 
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In a previous paper, we showed that the decay rate of a muon is only slightly affected by the 
presence of a circularly polarized laser and we gave an analytic expression for the correction. In this 
0^ , paper, we present the analytical result for the case of a linearly polarized laser. Again the effect of 

. the laser is small. 

o 
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INTRODUCTION 



Previously, some attempts have been made to find the change in the decay rate of the muon whenever a strong 
laser field is present. Liu, Li, and Berakdar [l| (LLB) tried to calculate this change for a strong linearly polarized 
laser, using an approximated electron wavefunction combined with numerical calculations. They found a large 
modification of the lifetime, as much as an order of magnitude. Narozhny and Fedotov challenged this result, 
arguing a small modification of the lifetime through a brief calculation 0, 0] . 

Recently, we showed the full analytical calculations of the decay rate of the muon in the presence of a strong 
circularly polarized laser field, finding only small (explicit) corrections to the unperturbed decay rate 

Although one might not intuitively expect major alterations of this conclusion once a laser with a different 
polarization is used, it proves to be worth going through the calculations for a strong linearly polarized laser. 
As we will show, this calculation is much more complicated and tedious. The electron wavefunction will involve 
two different exponents of the sine function, resulting in triple summations of the Bessel functions; in contrast 
with the circularly polarized case, where we had only one exponent of the sine function, and consequently, only 
one summation of the Bessel functions Q. Also, this calculation will settle any possible disagreement among the 
community regarding the difference in laser polarization, as well as providing a reference for the future work. 



THEORY 



X 

Muon decay is represented by 

M~0P) — » e~(p) + i/ e («i) + i^fe) , (1) 

where the arguments label the associated momenta. The wavefunction of the relativistic electron in a electro- 
magnetic field can be found by solving the Dirac equation with the electromagnetic potential term present. It 
was first presented by D. M. Volkov in 1935 0,0, 

(0 - e4 - m) tp e (x) =0. (2) 

For a linear polarization, is given by 

A»{x) = a^cos k ■ x (3) 

a" = (0,5/w) (4) 



* Electronic address: farzinni@msu.edu 

t Electronic address: dicus@physics.utexas.edu 

* Electronic address: repko@pa.msu.edu 

8 Electronic address: tinsley@hendrix.edu 



2 



with £ the amplitude of the electric field; k • £ = 0. We choose the photons to propagate along the z-axis, 
k^ = (ui, 0, 0, uj), and the Volkov solution gives 



„2 2 



il> e (x) = ( I + t^T ) u{p)e- iq - x + i ^ sin2k - x e- iS ^ sink - x (5) 
V 2p ■ k J 

with the electron's effective momentum and mass 

q " =P "-i^k k "> rn 2 =ml-^f (6) 

where too is the rest mass of the unaffected electron, mo = 0.511 MeV. Also, note that q ■ k — p ■ k. 

Following the standard S— matrix theory, the matrix element for the reaction, Mfi, may be extracted from 

S fi = -t JL J [Vv, 7 a (1 - 7s) Vv] $el X (1 - 75) 1> v .]#x = -^(2^) 4 <5 (4) (pf - Pi) Mfi . (7) 

In order to evaluate the left hand side of ([7|), the following generating function for Bessel functions from 9.1.41 
of 1 



e 2 



§z(t-l/t) 



n— — 00 




may be used to deduce the following (general) relations 

xe -i£ S m(z-M = e- mz x ^ i (J„+i(Oe l( " +1)0 ° + J n -i(0e <(n ~ 1) *') > ■ (9) 

»=-<» 1 57 - J„-i(Oe <(n - 15 *°) J 

Using (jnj) , and after some algebra, from |(7J) we find the matrix element for a given value of I and L (one index 
per exponent in |[SJ)) 

G 

= -= w(g 2 )7A (1 - 75) u{P)u{p) [A + Ai^] 7 A (1 - 75) w(?i) (10) 

where 

A = J L {B)J t [D) (11) 
A, = ^J L (B)[j e+1 (D) + J e ^(D)] (12) 

and 

D _ ep- a ^ B _ e 2 a 2 ^ 

p ■ k ' 8p ■ k 

The momentum conservation is then 

P» + {2L + ()k» = q^ + q? + q% . (14) 

Note that L has a prefactor 2, coming from the argument of sine in ([5]). Following the standard procedure of 
determining the decay rate, described for example in |9j, we see that once (O is squared, there are four summation 
indices, L, L' , £, and £'. One of the two delta functions, however, takes care of one of the summation indices 

t = 2(L-L')+£, (15) 
hence, we are left with three summation indices. The total decay rate is then given [j| by 

r= £ r L , LV (16) 

L,L'.e 
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where for each set of indices 




T w= III/ Y^5^(P- q - qi - q2 + (2L + e)k)^ Ml, 2{L _ L>)+i M L ,e. (17) 



s„,s,,s„ 



Using FORM [Tfj, all the traces in (jTTJ) can be evaluated, and we may perform the integration over the neutrino 
momenta by using the well-known relation 

/ ^i^i Si{Q - qi - q2)qC ^ = ^(QV^Q^MQ 2 ), (18) 
with Q = P - q + (2L + I) k. These then reduce ([17]) in muon rest-frame - to 

T W = J dE\ q \ J dzO(Q*)T L , L , e (19) 

where M is the muon rest mass, and E is q°. The object Tj, ^ £ is the square of the matrix element (|10[) . summed 
over the spin and integrated over the neutrino momenta, with relation (|15| implied (see APPENDIX for the 
explicit expression.) 

Turning off the laser corresponds to a = 0, which in turn, from (|13p . sets D and B equal to zero. In that 
case, given the properties of the Bessel function J„(0) = for and Jq(0) — 1, only the first line of the 
expression for Tl.l'.^ (see APPENDIX) is non-zero, and after integration, the unaltered decay rate of the muon 
is recovered 



r° = ^Sr (20) 



G 2 M 5 

1927T 

where, as in terms proportional to the electron mass have been neglected. 

As explained in our previous treatment Q , in (|19p the limits of integration are determined by the function, 
and the integration separates into two parts 

dE / dzQ{Q 2 ) = dE dz+ dE dz (21) 

J Jm J-l JM^ Jz L (E) 

with zl(E) coming from the the condition Q 2 > 

M 2 + 2M(2L + £)uj-2E[M+{2L + £)uj} 
ZL{E) = 2{2L + l)uE ' (22) 

Since the first term of |2T|) is independent of the summation indices, substitution of ([21]) into (fl9|) allows us 
to sum before integrating over this first term. In order to perform the summation per index of Ti,,u,lt we 
use X^-oc Jn{z)J n +k{z) = Jfe(0) (from 9.1.75 of @), the recurrsion relation for Bessel functions, n J n (z) = 
|(J n +l(z) + J„_i(z)), and J-„(z) — (—l) n J n (z). We note that (L — L') can only equal integer values; this 
notion forces many of the summation terms to be zero. After some tedious calculation of the summations, the 
integration may be performed, resulting in (neglecting the electron mass terms) 

r = r°{l + ^(-^+ln^]-^l(6-lQl^ + ^] }■ (23) 




As we will see, the term proportional to M 2 /m 2 cancels with a similar term coming from the second part of the 
integration in (|2"Tj) . 
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The second integration term can also be tackled in the same way described in [4|. The indices of the Bessel 
functions are limited by I < D, and L < B, as Bessel functions become very small once the index exceeds the 
argument. Therefore, given typical values such as ea ~ 10~ 4 MeV and ui ~ 1 eV, (2L + is always much less 
than 1 MeV. Hence, the corrections from this integration term are small, as the range of integration is small. 
The integration can, thus, be expanded in a Taylor series, and after keeping only the first non-zero term, the 
correction becomes [J| 

TC = T °J^ dz £ (2L + *)T W (24) 



where Tl l' e is the same as Tl.l',i but without the prefactor of 128 G 2 . Note that z L {^-) = -1 from ^ 7 and 
hence, once again we are allowed to switch the order of summation and integration. Although the calculation of 
(|24l) is lengthy and extremely tedious, it results in a quite simple expression 




Note the cancellation of the term proportional to M 2 /m 2 , once ([23)) and (|25|) are added. 



3. DISCUSSION 



Comparing the final result (|26p with that for the circularly polarized laser, 

n f e 2 a 2 r, uj u? s , M 2 3 lu u? 

Tcirc = r°U + 8 — 7T (— + 4—) In — + 2 — 

I M 2 l K M M 2 ' m 3 2M M 2 

there are some interesting differences. The corrections to the unperturbed decay rate are not the same in both 
polarization cases. There are terms proportional to e 4 a 4 /M 4 in a linearly polarized laser field and the correction 
in the linear case is larger than in the circular case because the \n(M/m) appears with a coefficient of 1 rather 
than a coefficient of u/M. Part of this difference is attributable to the fact that terms involving e a ^ 1& contribute 
to the circularly polarized case but do not contribute in the linearly polarized case. 

Our calculations in the present paper, as well as in the previous one show that the corrections to the 
undisturbed decay rate in both laser polarization cases are, however, consistently small, as one might expect 
intuitively. 



(27) 



4. APPENDIX 



The explicit expression for TljjjLi prior to the integration over the electron momentum and triple summation, 
and after simplification, is given by 

Tlx>, i =128G 2 x 

J J L (B)J L ,(B)J e (D)J 2(L _ L , )+i (D) \ (P - qf (P ■ q) + 2(P - q) ■ P (P - q) ■ q 



e 2 a 2 



+ ^[(P - q) 2 (P ■ k) + 2(P - q) ■ P (P - q) ■ k] 



(28) 



+ (2£ + 4L) J L (B)J L , (B)J e (D)J 2(L _ L , )+e (D) (P - q) ■ q (P ■ k) + (P - q) ■ k (P ■ q) 



+ (P-q)-P(q-k) 
+ (ML + 8L 2 + 2£ 2 ) J L (B)J L , (B)J e (D)J 2{L _ L>)+l (D) f (P • k) (q ■ k) 



eV 



2q ■ k 



(P-q)- k (P ■ k) 



JL(B)J L ,(B)J e (D)J 2{L _ L , )+e+1 (D)- 



[(q-a) (P-k)-(P-a) (q-k) 
+ i e aPlS k a q ^P s )] (Mi 2 - 4P ■ q + m ) 



+ (£ + 2L) J L (B)J L ,(B)J t (D)J 2(L _ u)+t+1 (D) 



2q ■ k 



[(q-a) (P-k)-(P-a) (q-k) 
+ i Sa^s^qPa-tp 5 )] (2P -k- q-k) 



+ J L {B)J L ,{B)J t {D)J 2{L _ L , )+i _ 1 {D) 



Aq ■ k 



+ (£ + 2L) J L {B) JuiB^D) J^l-^+i^D) 



bracket ([31 

e 



bracket 



J L (B)J L ,(B)J i+1 (D)J 2{L _ L , )+e (D)- 



2q ■ k 

[(q-a)(P-k)-(P- a) (q ■ k) 



(£ + 2L) J L (B)J L ,(B)J t+1 (D)J 2(L _ L , )+t (D) 



i e afhS k a q fi <i t P & )] (3M 2 - 4P ■ q + m 2 ) 
[(q-a) (P-k)-(P-a) (q-k) 
- i e aMS k a q l3 a < P s )\ (2P -k-q-k) 



2q ■ k 



+ Jl(B)J l ,(B)J 1 _ 1 (D)J 2{l _ l , )+1 (D) 



Aq ■ k 



bracket ([3"5]) 



+ (£ + 2L) J L (B) J L ,(B)J e ^(D) J 2[L _ L , )+t (D) 



2q ■ k 



bracket (j36j) 



Jl(B)J l ,(B)J 1+1 (D)J 2[l „ l , )+1+1 (D) 



e 2 a 2 



(P - q) 2 (P ■ k) + 2(P - q) ■ P (P - q) ■ k 



(I + 2L) J L (B)J L ,(B)J e+1 (D)J 2{L _ L , )+e+1 (D) 



e 2 a 2 
2q ■ k 



(P -q)-k(P-k) 



J L (B)J L ,(B)J l+1 (D)J 2[L _ L 



(P - qf (P ■ k) + 2(P - q) ■ P (P - q) ■ k 



(I + 2L) J L (B)Jv{B)Jt +1 {p)J 2{L - I/)+l _. x {D) 



2q ■ k 



(P-q)-k(P-k) 



J L (B)J L ,(B)J e . 1 (D)J 2{L _ L , )+e+1 (D)- t 



(P - qf (P ■ k) + 2(P - q) ■ P (P - q) ■ k 



(£ + 2L) J L (B)J L ,(B)J^ 1 (D)J 2(L ^ L , )+l+l (D) 



e 2 a 2 - 



2q ■ k 



(P-q)-k(P-k) 



J L (B)J u (B)J l _ 1 (D)J 2{L _ L , )+l _ 1 (D)- i 



■ ■ k 



(P-q) 2 (P-k) + 2(P-q)-P(P-q)-k 



(£ + 2L) J L (B)J L ,{E)Jt- i .(D)J 2{L _ IJ)+t _i(D) 



e 2 a 2 r 



2q ■ k 



(P-q)-k(P-k) 



(29) 
(30) 

(31) 

(32) 
(33) 
(34) 

(35) 

(36) 
(37) 

(38) 
(39) 
(40) 
(41) 
(42) 
(43) 
(44) 
(45) 
(46) 



where "bracket(*)" means the expression in the bracket of the corresponding "^''-equation, and £0123 = 1- 
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